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Forward-backward correlations of charged-particle multiplicities in symmetric bins in 
pseudorapidity are studied in order to gain insight into the underlying correlation structure 
of particle production in Au+Au collisions. The PHOBOS detector is used to measure 
integrated multiplicities in bins defined within |?7| < 3, centered at 1] and covering an 
interval Arj. The variance of a suitably defined forward-backward asymmetry variable 
is calculated as a function of rj, Arj, and centrality. It is found to be sensitive to short 
range correlations, and the concept of "clustering" is used to interpret comparisons to 
phenomenological models. 



1. Introduction 

It is often presumed that particle correlations are only "short-range" in rapidity. Short 
range correlations have been observed at all energies [Ej, in two-particle correlation mea- 
surements. However, they have never been shown to be the only source of correlation in 
multiparticle production. Indeed, single-particle distributions of inclusive charged parti- 
cles produced in heavy ion collisions reveal a distinct "trapezoidal" structure stretching 
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across the full rapidity range available in these reactions. More importantly, there is no 
evidence of an extended boost invariance, which would imply independent emission from 
different rapidity regions. Rather, two non-trivial effects are visible in the centrality de- 
pendence of particle production, which are apparently long-range in rapidity. Integrating 
over the full phase space reveals that particle production is linear with the number of par- 
ticipants This occurs despite a significant change in the shape of the pseudorapidity 
dependence as a function of centrality, which happens to be collision-energy indepen- 
dent in the forward region [Ej. All of this suggests that charged-particle production at 
mid-rapidity is highly correlated with particle production in the forward region, which 
naturally begs the question of the underlying structure of the single-particle distributions. 
In these proceedings, we discuss how to take first steps in this direction via the study of 
forward-backward multiplicity correlations. 

2. Forward-Backward Correlations from Fluctuations 

In this work, particle correlations are studied by the event-by-event comparison of the 
integrated multiplicity in a bin defined in the forward (77 > 0) region, called Np, centered 
at 7] with width At], with the multiplicity measured in an identical bin defined in the 
backward hemisphere (called 'W^"). With these definitions, one can construct the event- 
wise observable C = {Np — Nb)/\/Np + Nb, and measure its variance for a set of 
events with nominally similar characteristics (e.g. centrality). If particle sources tend 
to emit into the forward or backward region, such that the partitioning is binomial, this 
leads to 0"^ = 1, since a'^^Np — Nb) = Np + A''^. Short-range correlations arise if the 
objects emitted into either hemisphere tend to break into k particles, each of which stay 
close in rapidity (e.g. due to isotropic emission). If each cluster decays into exactly k 
particles, then C ^/kC, so cr^ — fccr^. Such intrinsic short-range correlations have in 
fact been measured in p + p and p + p experiments by direct construction of the 2-particle 
correlation function in 77 By varying Arj and measuring the ratio 4ap/{Np) as a 

function of Np + Nb, they found an "effective" cluster multiphcity (/ce// = {k) + o'l/ {k)) of 
approximately 2 charged particles, which exceeds the value of 1.5 estimated for a resonance 
gas [E]. 

3. Experimental Setup and Analysis Method 

The data analyzed here were taken with the PHOBOS detector [|6j during Runs 2 and 

4. in 2001 and 2004, respectively. The pseudorapidity acceptance was restricted to that 
of the "Octagon" detector, which is a tube of silicon sensors covering Ir^l < 3 and full 
azimuth except for a region near midrapidity. To simplify the analysis for different values 
of rj and Arj, only the regions of the detector with complete rapidity coverage were kept, 
restricting the total azimuthal acceptance to A0 = vr, in four 45-degree wedges [[7j. The 
multiplicity in each bin is estimated event-by-event by summing up the angle-corrected 
deposited energy of all hits and then dividing by the average energy per particle [IH]. 
Beyond the usual lower threshold to define a hit in the silicon, an ?7-dependent upper 
bound of the deposited energy per hit is applied to reduce the effect of slow secondaries 
on the fluctuations. The average (C), calculated as a function of rj, centrality, and event 
vertex, is subtracted event-by-event to correct for gaps in the Octagon (a process which, 
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according to simulations, leaves the fluctuations unaffected). 

To provide information that can be directly applied to models, a procedure was de- 
veloped to estimate and remove the detector effects from the raw measured value of 
{(Jcraw) by using Monte Carlo (MC) simulations of the PHOBOS apparatus. The basic 
idea is to assume that Ccraw — ~^ ^lev where a^^^ is the contribution from detector 
effects, and use the MC to subtract it on average, leaving behind only the physical corre- 
lations. However, it is found that there are several sources which contribute differently as 
a function of t] and combine in quadrature to a nearly-constant value over the pseudora- 
pidity range covered by the Octagon [|S]. These are corrected on average by calculating 
0"^ using a modified HIJING simulation with all intrinsic correlations destroyed, which 
allows a direct estimation of cr^g^. There remains a subdominant correlation of ajg^ with 
cr^ which is also removed. Finally, we correct for using only half-azimuth acceptance by 
cr^ — >■ 2(cT^ — 1) + 1, a formula obtained using MC simulations, assuming that the limited 
acceptance only affects the clustering properties of particle production. 

Systematic uncertainties were calculated by varying several assumptions about the es- 
timation of o"c< and found to be around Aa^ ~ 0.1. The bin-to-bin variation of the 
systematic error calculation was averaged over rj to reduce fluctuations in the error deter- 
mination procedure. 

4. Results 

After correcting for detector effects, the results on o"^ can be directly compared with 
model calculations based on charged primary particles. We have focused mainly on HI- 
JING [ ||9j and AMPT [ ^] , which have been used to describe various features of heavy 
ion collisions at RHIC energies. 

The first set of results concern the rj dependence, for forward and backward bins that 
are Arj = 0.5 units wide, as shown in Fig. ^ In this case, HIJING and AMPT already 
show some interesting differences. For peripheral (40-60%) events, both models have 
a similar magnitude and a monotonically-rising t] dependence. Central events show a 
substantial difference extending over most of the rapidity range, with AMPT showing a 
systematically smaller value of cr^'- This may be due to initially produced clusters being 
destroyed by the hadronic rescattering stage. Also, it is observed that both data and MC 
show cr^ ~ 1 at ?7 = 0. Although this is suggestive of the effect predicted by Jeon et al [ 
ITT] , it can also be explained if clusters produced at r] = emit particles into both the Np 
and Nb side, inducing an "intrinsic" long-range correlation that decreases a^. 

The next set of results is the Ai] dependence of a^, with a fixed rj = 2.0, as shown 
in Fig. n One sees that for both peripheral and central data, rises monotonically 
with increasing pseudorapidity interval. This can be explained in the context of cluster 
emission: by increasing At], one increases the probability of observing more than one 
particle emitted from a single cluster in either Np or Nb- Clearly, the rate of change with 
At] should reflect the full cluster distribution (both (k) and ak), but these studies have not 
yet been done. It is striking that the peripheral data has already reached cr^ ~ 3 for the 
largest At] in peripheral events, while this number is closer to 2 for central data. Finally, 
it is interesting that neither HIJING nor AMPT can explain both the centrality and Ar] 
dependence simultaneously, as HIJING reproduces the central data but has no centrality 
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Figure 1. (upper) for a fixed Arj = 0.5 as a function of t]. (lower) as a function of 
Arj for a fixed bin center at rj = 2.0. 



dependence, wliile AMPT is lower than the central data, but also sees a decrease in the 
overall magnitude moving from peripheral to central events. 



5. Conclusions 

In conclusion, measurements of forward-backward fluctuations provide insight into the 
structure of long and short range correlations in pseudorapidity space. The new PHOBOS 
data for 200 GeV Au+Au collisions are now fully corrected for detector and background 
effects, so direct comparisons can be made to phenomenological models. We see signif- 
icant short-range correlations at all centralities and pseudorapidities, instead of just at 
mid-rapidity. There is a non-trivial centrality and rapidity dependence of these corre- 
lations, in both 7] and A?^. Finally, neither HIJING nor AMPT reproduces even the 
qualitative features, but the way in which they fail to do so may well provide information 
on the underlying physics. In particular, more theoretical attention should be paid to the 
properties of "clusters" required to explain our data. Jeon et al have proposed, in Ref. [ 
ITT| . that the formation of a QGP near mid-rapidity should destroy any sort of cluster 
structure seen in p-|-p, and thus lead to a reduction in a^, beyond the detector-related 
effects discussed here. The data shown here should provide means to study such effects, 
or set upper limits on their occurrence. 
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